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ABSTRACT 

We present Atacama Large Millimeter/submillimeter Array (ALMA) 1.3 mm continuum images of the as¬ 
teroid 3 Juno obtained with an angular resolution of 0.042" (60 km at 1.97 AU). The data were obtained over a 
single 4.4 hr interval, which covers 60% of the 7.2 hr rotation period, approximately centered on local transit. 

A sequence of ten consecutive images reveals continuous changes in the asteroid’s profile and apparent shape, 
in good agreement with the sky projection of the three-dimensional model of the Database of Asteroid Models 
from Inversion Techniques. We measure a geometric mean diameter of 259±4 km, in good agreement with 
past estimates from a variety of techniques and wavelengths. Due to the viewing angle and inclination of the 
rotational pole, the southern hemisphere dominates all of the images. The median peak brightness temperature 
is 215± 13 K, while the median over the whole surface is 197 db 15 K. With the unprecedented resolution of 
ALMA, we find that the brightness temperature varies across the surface with higher values correlated to the 
subsolar point and afternoon areas, and lower values beyond the evening terminator. The dominance of the 
subsolar point is accentuated in the final four images, suggesting a reduction in the thermal inertia of the re- 
golith at the corresponding longitudes, which are possibly correlated to the location of the putative large impact 
crater. These results demonstrate ALMA’s potential to resolve thermal emission from the surface of main belt 
asteroids, and to measure accurately their position, geometric shape, rotational period, and soil characteristics. 

Subject headings: minor planets, asteroids: general — minor planets, asteroids: individual (3 Juno) — planets 
and satellites: surfaces — techniques: interferometric 
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1. INTRODUCTION 

Discovered in 1804, Juno was the third main-belt asteroid 
identified, following Ceres and Pallas. The first reasonably 
accurate measurement of Juno’s diameter was performed with 
filar micrometers on the Great Lick Refractor (Barnard 1895) 
and the Yerkes 40-inch Refractor (Barnard 1900), yielding 
a value of 193 ±20 km (see also Dollfus 1971). A mod¬ 
ern measurement of its physical cross section came from 15- 
station stellar occultation data, yielding a mean diameter of 
267±5 km with a significant ellipticity (Millis et al. 1981). 
Optical speckle interferometry soon produced a size measure¬ 
ment consistent with the occultation result (Baier & Weigelt 
1983). Like most asteroids, Juno’s light curve is double- 
peaked with two maxima and two minima (e.g., Birch & 
Taylor 1989), indicative of a non-spherical shape. Based on 
light curve inversion, Juno has a unique rotational pole that 
is significantly tilted with respect to the ecliptic (Magnus- 
son 1986; Dotto et al. 1995), and its period of 7.209531 hr 
is known to high accuracy (Kaasalainen et al. 2002). This in¬ 
formation, combined with recent, near-infrared adaptive op¬ 
tics (AO) imaging led to a triaxial ellipsoid model with axis 
lengths of 297, 233, and 222 km (Drummond & Christou 
2008). A three-dimensional model with 2036 faces and 1020 
vertices based on a combination of the historical optical light 
curves and two occultations (Durech et al. 2011; Kaasalainen 
et al. 2002) is hosted by the Database of Asteroid Models from 
Inversion Techniques (DAMIT; Durech et al. 2010). 

Juno is a member of the S-class of asteroids (Chapman et 
al. 1975), which have a stony composition of iron-bearing sil¬ 
icates and metallic iron as inferred primarily from their 1 pm 
spectral absorption feature (e.g., Gaffey et al. 1993a, and ref¬ 
erences therein). The optical Small Main-Belt Asteroid Spec¬ 
troscopic Survey (SMASSII) assigns it subclass Sk as a tran¬ 
sition object toward the K-class, which exhibits a shallower 
1 pm feature (Bus & Binzel 2002). Infrared Space Observa¬ 
tory spectra of Juno show an 8-11.5 pm feature that is consis¬ 
tent with the laboratory measurements of the silicate minerals 
pyroxine and olivine (Dotto et al. 2000). Evidence for sur¬ 
face features on Juno have been suggested by the variation 
as a function of rotation angle of its optical colors (Degewij 
et al. 1979; Schroll et al. 1981) and linear polarization (Shi- 
nokawa et al. 2002; Takahashi et al. 2009), and by a sequence 
of optical AO images, which suggested a large impact crater 
(Baliunas et al. 2003). Somewhat surprisingly, there are no 
published images of Juno from the Hubble Space Telescope 
(Dotto et al. 2002), and there have been no spacecraft encoun¬ 
ters as yet. 

As a powerful new tool in the study of Solar System bodies, 
the Atacama Large Millimeter/submillimeter Array (ALMA; 
Hills et al. 2010) will be able to map the shape and surface 
temperature distributions of hundreds of main belt asteroids 
and Jupiter Trojans (Busch 2009; Lovell 2008). The reason is 
twofold. First, the absorption length of (sub)millimeter pho¬ 
tons (Campbell & Ulrichs 1969) on asteroid surfaces is com¬ 
parable to the thermal skin depth of the diurnal wave (typ¬ 
ically a few to 10 wavelengths; Spencer et al. 1989); thus, 
these wavelengths are well-matched to probe the thermal re¬ 
sponse of this material and should provide information on the 
thickness, structure and nature of the regolith (Chamberlain 
et al. 2007; Lagerros 1996). Indeed, the recent flybys of the 
main belt asteroids 21 Lutetia and 2867 Steins by the Rosetta 
mission have yielded important measurements of the ther¬ 
mal inertia and emissivity of their surface material (§4) us¬ 


ing its (sub)millimeter radiometers (Gulkis et al. 2010, 2012). 
Second, ALMA has exquisite continuum brightness tempera¬ 
ture sensitivity at small angular scales. For example, the 50- 
antenna, full-bandwidth sensitivity in one minute at 300 GHz 
(1 mm) at the highest angular resolution (13 mas) is 10 K. 
Since their physical temperatures are typically 100-200 K, 
ALMA can effectively image these bodies down to linear res¬ 
olutions of 10 km (at a distance A = 1 AU) with high signal- 
to-noise ratio (S/N), enabling the use of the powerful tool of 
self-calibration. As a significant step toward demonstrating 
this capability, in this Letter we present the first millimeter 
wavelength images to resolve the surface of Juno, which were 
obtained during the recent campaign to commission ALMA’s 
long baseline capabilities (ALMA partnership et al. 2015). 
The availability of the DAMIT model provides an excellent 
test of ALMA’s imaging performance, while the resulting im¬ 
ages provide new details on the surface conditions of Juno. 

2. OBSERVATIONS 

An approximately 53-minute length scheduling block (SB) 
to observe the 1.3 mm (233 GHz) continuum emission from 
Juno was executed five consecutive times on 2014 Oct 19 
starting at 09:15 UT (43 min before local sunrise) and end¬ 
ing at 13:38 UT. Four spectral windows were used, each with 
2 GHz bandwidth, 128 channels and dual polarization. Center 
frequencies were 224, 226, 240, and 242 GHz. All necessary 
calibration observations were performed in each execution of 
the SB. An additional focus measurement and adjustment was 
performed prior to the fifth execution (two hours after sun¬ 
rise) as per normal operations. The SB included an exter¬ 
nal ephemeris with 4 min sampling obtained from Jet Propul¬ 
sion Lab (JPL) Horizons 17 , which reports a 3 a uncertainty of 
60 mas in right ascension and 26 mas in declination. 

Calibration and imaging was performed in CASA 18 version 
4.2.2. The complex gain calibration cycle time was 68 s, 
with Juno being observed for 48 s and the gain calibrator 
(J0757+0956) for 15 s (5.7° away). Data from 27 to 31 an¬ 
tennas were used, ranging in projected baseline length from 
0.02 to 10 MA (26 m to 13 km). The zenith precipitable water 
vapor varied from 1.4 to 1.6 mm. Bandpass and flux cali¬ 
bration is based on observations of the quasar J0750+1231 
in each SB. This quasar is an ALMA calibration grid source 
(van Kempen et al. 2014) for which a linear interpolation in 
frequency from the measurements nearest in time at 3 mm (7 
days) and 0.87 mm (18 days) yields an assumed flux density 
of 0.64±0.04 Jy at 233 GHz with a spectral index of -0.66. 
The calibrator flux density measurements were stable over 
many weeks and we estimate our flux scale to be accurate 
to 6%. The mean flux density derived for the gain calibrator 
was quite consistent across the five executions, differing by a 
maximum of 2.4% from the weighted mean of all executions 
(0.5916 ± 0.0007 Jy). Thus, there was very little decorrela¬ 
tion on timescales shorter than the integration time (1.92 s). 
As shown in Table 1, the span of observations was approxi¬ 
mately centered on the time of transit of Juno. 

Following calibration, the uv data from each execution was 
split into two halves, with the duration and time on source 
of each half being ^18 min and ~10 min, respectively. The 
ten resulting datasets (epochs) were imaged individually us¬ 
ing a robust weighting parameter of zero. Phase-only self¬ 
calibration was then performed, initially with a time interval 

17 http://ssd.jpl.nasa.gov/horizons.cgi 

18 http://casa.nrao.edu 
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of 300 s, followed by a refinement with a time interval of 15 s. 
The final execution showed phase-cal solutions of somewhat 
larger magnitude and higher variability than the rest, with a 
loss of some antennas on the outermost pads. Amplitude self¬ 
calibration on a timescale of 300 s was then performed. The 
final images were constructed using multi-scale clean (Rau 
& Cornwell 2011) with deconvolution scales of 0, 5, and 15 
pixels to avoid image artifacts caused by the clean instabil¬ 
ity that can occur when modeling an extended source using 
only delta functions. The image dynamic range after self¬ 
calibration was 120. This improvement factor of 2-6 indi¬ 
cates that a significant residual phase error was present af¬ 
ter normal calibration. Thus, it is important to realize that 
any usage of calibrated visibilities for direct modeling (e.g., 
Viikinkoski et al. 2015) must either apply the imaging self¬ 
calibration solutions or include antenna-based phase solutions 
as a model parameter to be solved (Hezaveh et al. 2013). As 
expected when the initial self-calibration model has high S/N, 
the intensity-weighted centroid of the source before and after 
self-calibration is consistent to within a fraction of the (5 mas) 
pixel size (< 0.2 in most images, and < 0.3 in the final two 
images). Because we used an accurate VLBI position for 
the gain calibrator (07:57:06.64296, +09:56:34.8525; Lanyi 
et al. 2010), we expect our images to follow ALMA’s mea¬ 
sured astrometric performance (3 mas; ALMA partnership et 
al. 2015) 19 . The images vary in the size of the synthesized 
beam from 31.8 x 23.7 mas to 41.8 x 36.1 mas, with a mean 
position angle of 30° ±11°. These images are publicly avail¬ 
able from the ALMA Science Verification page 20 . To present 
matched resolution images, we smoothed the images with a 
two dimensional elliptical Gaussian to obtain a circular beam 
of 42.0 mas, which at the Earth-Juno distance of A = 1.97 AU 
corresponds to 60 km. No near-field correction was applied 
to the uv data, but Juno was beyond the Fraunhofer distance 
for the longest baseline, and the self-calibration process may 
mitigate the residual effects. 

3. RESULTS 
3.1. Viewing geometry 

Juno’s orbit has a mean radius of 2.67 AU with signif¬ 
icant eccentricity (0.25) and inclination (13°). As shown 
in Figure 1, the true anomaly was 37.5° during our obser¬ 
vations, with a solar phase angle of -28° yielding an illu¬ 
mination percentage of 94%. We have plotted the orienta¬ 
tion of the rotational pole inferred from parametric blind de- 
convolution of near-infrared AO images (ecliptic longitude 
A = 118° and latitude f3 = +30°) which has an uncertainty of 
13° (Drummond & Christou 2008). The direction of Juno was 
A = 124.7°, (3 = -13.2°, thus our viewing angle was 48° from 
the polar axis but nearly coplanar, with the southern hemi¬ 
sphere dominating the view, as shown in the inset of Figure 1. 
The angular smearing due to the 15° of axial rotation during 
each 18 min image leads to 21 mas of linear smearing at the 
mean radius, which is half the beam size, leading to < 12% 
loss in resolution. The mean epochs of the images are listed in 
Table 1, along with the corresponding rotational phases com- 

19 The ALMA control system does not account for the finite distance of 
solar system targets when computing the gravitational deflection of light by 
the Sun to apply to the astrometric ephemeris, which causes the phase transfer 
from the quasar to introduce a systematic position error of « 1.4 mas for these 
observations of Juno according to the formulas of Cowling (1984). 

20 http://almascience.org/alma-data/science-verification. Additional de¬ 
tails including the calibration and imaging scripts are available from 
http://casaguides.nrao. edu/index.php?title=ALMA2014_LBC_SVDATA. 


TABLE 1 

Parameters of of the 1.3 mm ALMA images of Juno 


Image 

# 

Epoch 

(MJD) 

Elapsed time 
(minutes) 

Elevation 

(°) 

Rotation phase 

(fiV 4>2 b 

0 

56949.39167 

0 

53.2 

0.33 

0.18 

1 

56949.40417 

18 

55.7 

0.37 

0.22 

2 

56949.42917 

36 

59.5 

0.45 

0.30 

3 

56949.44167 

73 

60.5 

0.49 

0.35 

4 

56949.46667 

110 

60.5 

0.58 

0.43 

5 

56949.48056 

128 

59.4 

0.62 

0.47 

6 

56949.51111 

173 

54.7 

0.72 

0.58 

7 

56949.52500 

192 

51.9 

0.77 

0.62 

8 

56949.54861 

226 

46.1 

0.85 

0.70 

9 

56949.56111 

244 

42.6 

0.89 

0.74 


a Phase with respect to t/’o of the Drummond & Christou (2008) triax- 
ial model using the rotation period of 0.300397125 days (Kaasalainen et 
al. 2002). In the 8.5-year interval, the phase accumulated uncertainty is 
± 0 . 01 . 

b Phase with respect to AO observations of Baliunas et al. (2003) (zero 
point taken to be 50371.24167) using the same rotation period. In the 
18-year interval, the accumulated phase uncertainty is ±0.02. 

puted with respect to the zero time point of the Drummond & 
Christou (2008) triaxial model, and with respect to the optical 
AO observations of Baliunas et al. (2003). 

3.2. Images 

The ten images of Juno are shown in Figure 2. The absolute 
position was measured in each image by computing the cen¬ 
troid of all pixels above the 5 a level, where cr, the image rms, 
was defined by an annulus surrounding the object (see Table 2 
notes). These pixels were weighted uniformly to avoid influ¬ 
ence of surface brightness variations across the face of the 
object. The difference between the centroid position and the 
image phase center yields the observed offset from the JPL 
ephemeris, which is stable in right ascension at ~ +60 mas, 
but slowly varying in declination. The integrated flux den¬ 
sity of the source was measured by integrating over all pixels 
above the 3 a level. The peak positions, intensities, and the 
corresponding Planck brightness temperatures (T B ) are also 
listed in Table 2. Note that these 7 b differ from the Rayleigh- 
Jeans approximation by +5.6 K. An estimate of the median 
brightness temperature across the surface was computed for 
each image by finding the median pixel intensity in the clean 
component model image and dividing by the solid angle of a 
pixel. The result is typically 10-15 K below the peak T B . 

3.3. Size and shape measurement 

In order to obtain a size measurement independent of the 
existing shape models, we consider the underlying source to 
be an elliptical disk with uniform brightness. First, we first fit 
the observed image from each epoch with a two-dimensional 
elliptical Gaussian, recording the position angle and major 
and minor axes as the target parameters. We then create a 
disk model image using the target parameters, but with the 
major and minor axes increased by 30% to account for the 
bulk of the size underestimate resulting from the Gaussian 
brightness profile. We then convolve this disk model image 
with a 42 mas beam to match the observations. Next we it¬ 
eratively refined the disk model parameters until an elliptical 
Gaussian fit to the convolved disk model image matched the 
target parameters. The resulting uniform brightness elliptical 
disk model parameters for each epoch are listed in Table 2. 
The geometric mean of the median of the major and minor 
axes of the ten elliptical disk models (0.181" ± 0.003") cor- 
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Fig. 1. — As viewed from the ecliptic pole, this diagram shows the alignment of Earth and Juno in their respective orbits on the three dates of observational data 
discussed: ALMA (this paper), Mt. Wilson Observatory (Baliunas et al. 2003), and Lick Observatory (Drummond & Christou 2008). The scale of heliocentric 
ecliptic coordinates (A, f$) is indicated by the dashed line circle. The three angles inside the dotted lines correspond to the sun-target-earth angle which describes 
the solar illumination phase. The nominal rotational pole (toward A = 118°,/5 = +30°; Drummond & Christou 2008) and its corresponding equator are drawn 
onto a spherical representation of Juno for reference. Juno’s perihelion at A = 58.3° is marked by the dot labeled “P”. The portion of the orbit above the ecliptic 
is shown by the thick line, which begins at the ascending node (A = +169.9°). Inset: The inset shows the point of view of the ALMA observations (in the ecliptic 
coordinate frame), which is dominated by the southern hemisphere. Juno’s south pole and equator are marked, as is the evening terminator and unlit side which 
reflects Juno’s solar phase angle and heliocentric latitude (/3 = -13.2° viewed from Earth) at the time of observation. 


TABLE 2 

Derived properties from the 1.3 mm ALMA observations 


Image 

# 

Centroid of Juno a 

Absolute position 
a (J2000) 6 (J2000) 

Offset 1, 

Aa, A<5 
mas, mas 

Integrated 
Flux density c,d 
mJy 

Peak 

Intensity d,e 
mJy beam -1 

Peak d 

K 

r B 

Median d,f 

K 

7b Peak 
Position 

mas, mas 

Elliptical disk model 
Parameters 8 
mas x mas (°) 

0 

08 14 46.8474 

+06 19 08.554 

+59, -5 

199.7±1.5 

17.0±0.09 

2223=1 

213 

-9, -5 

200x157 (+47±3) 

1 

08 14 47.8979 

+06 19 02.586 

+55, +4 

198.1±1.5 

16.8±0.09 

220±1 

208 

-20, +21 

199x159 (+37±3) 

2 

08 14 49.9975 

+06 18 50.632 

+52, +14 

199.23=1.5 

16.5±0.09 

216±1 

202 

-22, +36 

196x165 (+16±4) 

3 

08 14 51.0474 

+06 18 44.645 

+62, +13 

200.6±1.5 

16.6±0.09 

218±1 

200 

+8, -3 

192x169 (+0±6) 

4 

08 14 53.1449 

+06 18 32.674 

+63, +23 

202.9±1.5 

16.3±0.09 

2143=1 

197 

+7, +7 

191x175 (-47±8) 

5 

08 14 54.3100 

+06 18 26.017 

+63, +26 

201.8±1.4 

16.2±0.09 

213±1 

196 

+27, +34 

194x173 (-67±6) 

6 

08 14 56.8740 

+06 18 11.357 

+62, +31 

195.8±1.5 

15.7±0.09 

207±1 

186 

+38, +9 

200x170 (-99±4) 

7 

08 14 58.0396 

+06 18 04.685 

+55, +30 

196.1 ±1.5 

15.8±0.09 

208±1 

188 

+40, +10 

200x169 (-108±4) 

8 

08 15 00.0242 

+06 17 53.355 

+62, +47 

188.43=1.8 

16.3±0.11 

215±2 

187 

+33, +13 

192x167 (-126±5) 

9 

08 15 01.0759 

+06 17 47.339 

+64, +43 

181.6±2.0 

16.4=L0.12 

215±2 

179 

+31, +2 

187x168 (-137±6) 

Median 11 



+60.5, +24.5 

198.7±2.7 

16.4=L0.2 

215±3 

197 ±9 


194 x 169 


a See section § 3.2 for details of the centroid calculation. The systematic uncertainty is estimated to be 3 mas. All coordinates are in the ICRL of 
the ICRS (Epoch J2000.0). 

b Offset of centroid with respect to the phase center of the image, which follows the JPL Horizons ephemeris. 

c Flux density integrated over pixels in the region with intensity >3 cr, where cr (the image rms) is calculated in an annulus of radii 0.325" to 
0.575". The quoted uncertainty is: (Number of independent beams in the region) 0 5 x 3 cr, where cr is the image rms. 
d The systematic flux calibration uncertainty of 6% is not included. 
e The quoted uncertainty is the image rms. 

f The median brightness temperature (7 b) is computed using the Planck equation with the median pixel intensity in the clean model (within the 
central 0.15") and the solid angle per pixel. 

g Major axis, minor axes and position angle east of north. The fit uncertainty on the major and minor axes is 3 mas. 
h Where listed, the uncertainty is the median absolute deviation from the median. 


responds to d mm = 259 ± 4 km. This mean diameter is consis¬ 
tent with the size derived from 250 GHz single dish flux den¬ 
sity measurements under the assumption of unity emissivity 
(e v - 1; d mm = 253.4 ± 7.4; Altenhoff et al. 1994). The mean 
diameter is also in reasonable agreement with the triaxial geo¬ 
metric mean (250 ± 5 km) of the three axes of the Drummond 
& Christou (2008) triaxial ellipsoid model, the equivalent di¬ 
ameter of the Durech et al. (2011) model (252 =b 29 km), and 
the effective diameter measured by radar (265 ±30 km; Magri 


et al. 2007). The IRAS Minor Planet Survey (IMPS) inferred a 
radiometric mean diameter of 234 ±11 km, but IMPS diame¬ 
ters are systematically low compared to occultation diameters 
(Tedesco et al. 2004). We note that our assumption of uni¬ 
form brightness is somewhat unphysical, as the daytime sur¬ 
face temperature will vary as a function of local hour and lat¬ 
itude in more realistic thermal models (Lebofsky & Spencer 
1989). We estimate that our method will yield sizes that are 
^2-4% smaller than the mean physical size. Thus, further de- 
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Fig. 2.— In the two pairs of rows, the lower set of panels (with tick marks) are ALMA 1.3 mm continuum images of Juno in equatorial offset coordinates 
referenced to the phase tracking center corresponding to the JPL ephemeris, but with the centroid shifted by the values in columns 4 and 5 in Table 2. In each 
frame, the cross marks the position of the peak intensity. The position angle of the Sun was +101.5° east of north in equatorial coordinates, as indicated by the 
arrow. The beam size is indicated by the circle in the lower right-hand corner. The 3 a position uncertainty of the JPL ephemeris in each axis is indicated by 
the cross in the lower left-hand corner. The rotational phase values (<±) are from Table 1. The faceted images are DAMIT sky projection models (with artificial 
lighting) computed for the corresponding mean epoch of each image and are shown on the same angular scale (the models have a 12% uncertainty in scale). 


tailed comparisons of the millimeter images to any specific 
shape model should be performed in the context of a thermal 
model. 

3.4. Comparison to DAMIT shape model 

The model images shown in Figure 2 were obtained using 
the DAMIT online tool 21 (see also the Interactive Service 
for Asteroid Models 22 (ISAM); Marciniak et al. 2012). This 
tool provides a prediction of the projected appearance of the 
asteroid on the sky as viewed from Earth in equatorial co¬ 
ordinate orientation at any observed Julian date. The light 
travel time effect is taken into account (in this case 16.4 min¬ 
utes). The DAMIT prediction for the ESO 1.5 m speckle ob¬ 
servation of Baier & Weigelt (1983) is in excellent agreement 
with the shape of the observed image, obtained over 34 years 
ago. Likewise, the shape of the DAMIT prediction is in good 
agreement with most of the ALMA images. The images we 
show here use artificial lighting in order to show the full geo¬ 
metrical extent of the body, which will emit millimeter emis¬ 
sion that is only mildly modulated by solar illumination. In 
a qualitative sense, it is perhaps images 4 and 5 that are most 

21 http://astro.troja.mff.cuni.cz/projects/asteroids3D/web.php 

22 http://isam.astro.amu.edu.pl 


discrepant from the model in terms of ellipticity and orienta¬ 
tion. As for angular scale, the quoted accuracy of the DAMIT 
images is ±12%. The ALMA-derived major axes (Table 2) 
are systematically ~ 6% smaller than the maximum extent in 
the DAMIT images, however much of this difference could be 
due to the simplistic model used in § 3.3. 

3.5. Surface features 

To accentuate the variation of brightness across the object, 
we have created a model image of uniform brightness cor¬ 
responding to each of the ten images. We begin with the 
clean component model image, and compute the median value 
within the central 0.15" diameter (§ 3.2). We then place this 
value into all pixels of the clean model that are inside the 
40% level in the clean image, and place zero elsewhere. This 
approach leads to a model image comparable to the angular 
dimensions of Juno. We then convolve this model with the 
42 mas beam and subtract it from the clean image. This sub¬ 
traction enables the identification of areas of lower vs. higher 
brightness temperature, regardless of the relative calibration 
accuracy between the different images. The results are shown 
in Figure 3. In most of the frames there is a consistent pat¬ 
tern of the northwest edge being the coolest portion, which 
matches the location of the evening terminator. Also, the 
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warmest part of the image appears correlated with the subso¬ 
lar point. In the first five images, the afternoon area following 
the subsolar point is the warmest point. Meanwhile, in the 
last four images, the warmest point is at (or very close to) the 
subsolar point. The temperature contrast becomes particularly 
pronounced in the final two images as the surface median T B 
declines (Table 2). 

4. DISCUSSION 

In the simple equilibrium model, the expected disk- 
averaged brightness temperature, 7#(z/), of an asteroid is de¬ 
termined by its spectral emissivity, e v , and its mean equilib¬ 
rium temperature, T eq : 


Tb(v) = e u T eq 

(1) 

1 

II 

(2) 


where A is its bolometric Bond albedo, r is its heliocentric dis¬ 
tance in AU, and / = 329 K for non-rotating objects and 277 K 
for fast-rotating objects (Cremonese et al. 2002; Kellermann 
1966). The bolometric Bond albedo is a product of the bolo¬ 
metric geometric albedo ( p ) and the bolometric phase integral 
(< q ). In principle, both of these components represent integrals 
of wavelength-dependent quantities which must be measured 
and weighted by the solar flux spectrum (Hansen 1977), but 
they are often approximated as being wavelength indepen¬ 
dent. Measurements of p for Juno at optical wavelengths (py) 
range from « 0.13 (Hansen 1977; Brown et al. 1982) to 0.15 
(Morrison 1977; Zellner et al. 1977), while a value of ~ 0.22 
has been found at mid-infrared wavelengths (Ryan & Wood¬ 
ward 2010; Tedesco et al. 2004). The observed variation in 
optical brightness versus phase angle yields a slope parameter 
G = 0.17 dz 0.03 (Lagerkvist et al. 1992), which in turn yields 
qy = 0.41 ±0.02 from the relation of Bowell et al. (1989). 
Combining qy with the more recent measurements of pir, we 
will proceed with A = 0.09, which matches the result others 
have obtained by using the py values along with qy = 0.6 ap¬ 
propriate for the Moon. In any case, the dependence of T eq on 
A is quite mild. Assuming e v = 1, the equilibrium model pre¬ 
diction for T b for our Juno observations (r = 2.072 AU) is then 
188 K to 223 K, depending on /. Juno is a fast rotator, but its 
polar axis was pointed significantly in line with the Sun dur¬ 
ing the ALMA observations (Fig. 1). Therefore, an / value in 
the lower half of the range is likely to be appropriate, which 
is consistent with our measured median 7g bs (233 GHz) of 
197 ± 15 K. A previous single-dish estimate of Tg bs (345 GHz) 
was 146 ±48 K at r = 3.132 AU (Chamberlain et al. 2009), 
which scales to 180 ± 59 K at our smaller value of r. 

Moving beyond the equilibrium model, we next consider 
the Standard Thermal Model for asteroids (STM; Lebofsky & 
Spencer 1989). Because ALMA resolves the surface of Juno, 
we can compare the observed peak brightness temperature 
with the expected subsolar surface temperature (T ss ) in the 
STM. For Juno’s values of A and r, and its beaming factor of 
rj = 0.76 (Spencer et al. 1989), the expected T ss is 286 K. This 
value is significantly higher than the peak T B of 222 ± 13 K 
observed by ALMA, in contrast to previous mid-infrared mea¬ 
surements of Tss on Juno which are consistent with the STM 
prediction (Lim et al. 2005). If interpreted in the context of a 
model where all of the emission arises from the surface, such a 
discrepancy could be interpreted as an “effective” e u ~ 0.8 at 
1.3 mm. Effective emissivity is a quantity which can encom¬ 
pass many effects besides the physical emissivity of the ma¬ 
terial, including sub-surface sounding of deeper colder layers 


(Fornasier et al. 2013; Redman et al. 1998). Indeed, in terms 
of a physical model with temperature gradients in the sub¬ 
surface material of up to 50 K mm -1 , such a large temperature 
discrepancy is expected to be seen in millimeter wavelength 
emission, which arises from material at a range of depths, 
even when the bulk material has e u ~ 1 (see, e.g., Keihm et 
al. 2013). 

To interpret the enhanced brightness temperature of the sub¬ 
solar point in most of the Juno images, we consider heat con¬ 
duction in the regolith following the equations in Lagerros 
(1996). The thermal inertia (T) is a function of the surface ma¬ 
terial density (p), thermal conductivity (k), and specific heat 
capacity (c p ) of the soil, while the thermal skin depth (4) also 
depends on the angular rotation rate (cj), 

r = ^/Kpc p (3) 



Using the surface density of Vesta derived from its radar 
albedo (1.75 g cm -3 ; Chamberlain et al. 2009), along with 
values for k and c p in the porous lunar surface layer (2 x 10 -5 
W cm -1 K -1 and 0.6 J g -1 K -1 , respectively; Keihm 1984), 
yields T = 46 J m -2 s -0 5 K -1 and / s = 2.8 mm. Since 4 is 
only 2.2 wavelengths, we can expect the observed continuum 
emission to arise from a mix of solar-heated surface material 
and deeper unheated material. Thus, the correlation of Juno’s 
brightness temperature with the subsolar point is not surpris¬ 
ing. However, the fact that the brightest point moves from the 
subsolar afternoon area in the first half of the images to near 
coincidence with the subsolar point in the latter half of the im¬ 
ages, suggests a change in the soil properties with longitude. 
For example, in the case that k is constant with depth and lon¬ 
gitude, 4 would scale inversely with F due to changes in either 
p or c p . Thus, if soil with relatively lower values of p or c p ex¬ 
ists at these longitudes (0 \ ~ 0.7-0.9), it could have a lower 
thermal inertia, deeper skin depth, and consequently a greater 
proportion of the millimeter emission arising from heated ma¬ 
terial. A lower inertia across most of this side of Juno might 
also explain the lower median T B observed, as cooling would 
proceed more rapidly as the angle from the subsolar point in¬ 
creases. In any case, variations in thermal inertia across the 
surface of an asteroid are quite plausible, particularly in light 
of the detailed variations on Vesta reported from the Visual 
and Infrared mapping spectrometer on the Dawn spacecraft 
(Capria et al. 2014). 

Recently, two asteroids have been measured at 0.53 and 
1.6 mm by the Microwave Instrument for the Rosetta Or- 
biter (MIRO) during close encounters by the European Space 
Agency Rosetta spacecraft: the small (~6 km) object Steins 
(Gulkis et al. 2010), and the larger (^100 km) object Lute- 
tia (Gulkis et al. 2012). In contrast to Steins, which has a 
high, rock-like inertia (T = 450-850 J m -2 s -0 5 K -1 ) and 
e u = 0.85-0.9, Lutetia exhibits a very low thermal inertia 
(T = 20 J m -2 s -0 ’ 5 K -1 ) in the upper 1-3 cm much like the 
fine dust of the lunar regolith, with an emissivity consistent 
with reflection from a surface with dielectric constant of 2.3 
( e v = 0.958). These latter properties may be similar to what 
ALMA has seen on Juno, particularly in the latter half of the 
images. This result is perhaps not surprising in the context 
of the subsequent thermophysical modeling of Keihm et al. 
(2013), which finds that low thermal inertias and e v ~ 1 can 
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Fig. 3.— Residual images of Planck brightness temperature (7b) in equatorial offset coordinates created by removing a model image of uniform brightness 
from the images in Figure 2. The open circle indicates the location of the subsolar point. The color scale range is ±30 K with respect to the median values of 7 b 
in column 8 of Table 2. The south pole drawn is the same as in Figure 1, and the sense of rotation is clockwise. The different position angle here is due to the 
combination of the different coordinate system, the foreshortening effect of the heavily inclined pole, and the significant heliocentric latitude. 


fit the infrared to centimeter spectral energy distributions of 
the asteroids Ceres, Vesta, Pallas, and Hygiea. 

The possibility of a recent impact on Juno was raised by the 
detection of a region of reduced 934 nm brightness in AO im¬ 
ages correlated with a spatial “bite” feature on the limb (Bal- 
iunas et al. 2003). Unfortunately, the north/south orientation 
of the Baliunas et al. (2003) images is not specified. How¬ 
ever, judging from the DAMIT model images at that epoch, 
the AO images would appear to be oriented with south up be¬ 
cause in this case, a depression in the model images would 
map closely to the proposed crater near the limb of the fifth 
AO image. If so, then the crater is located near the north pole. 
As shown in Figure 1, the ALMA viewing angle of Juno dif¬ 
fers by 114° from the Baliunas et al. (2003) AO images such 
that the north pole is not visible. On the other hand, if the AO 
image is oriented with north up, as it is in similarly-acquired 
images of Vesta by a subset of these authors (Shelton et al. 
1997), then the crater would lie between the equator and the 
south pole, placing it at a latitude that crosses near the cen¬ 
ter of the ALMA view. The rotational phase of the fifth AO 
image (02 = 0.29) is close to that of ALMA image 3, thus 
a feature on the limb in the AO image would cross the sub- 
Earth point in the ALMA image 114°+90 o =204° later (i.e. 
at 02 = 0.86). This phase corresponds to ALMA images 8 
and 9. Those images do show the highest temperature con¬ 
trast with respect to the subsolar point, which could be con¬ 
sistent with a lower thermal inertia in the excavated material 
surrounding the crater. Clearly, future ALMA observations of 
complete rotations of Juno, preferably at multiple phases and 
wavelengths, will be necessary to explore this phenomenon 
further and develop an accurate, full-surface thermophysical 
model. 

5. CONCLUSIONS 

Our ALMA long-baseline observations of Juno provide the 
first ground-based images that significantly resolve the sur¬ 


face of an asteroid at millimeter wavelengths. They provide an 
independent set of size and shape measurements which con¬ 
firm our current knowledge expressed by the DAMIT and tri- 
axial ellipsoid models. Future ALMA observations of main 
belt asteroids, including both spatially-unresolved photomet¬ 
ric lightcurves (e.g., Moullet et al. 2010) and resolved images, 
can be used to test and refine the existing three-dimensional 
models. We note that ALMA can potentially achieve signifi¬ 
cantly higher physical resolution on Juno that these initial ob¬ 
servations offer. For example, a factor of three improvement 
(to 20 km resolution) would be possible by observing in the 
345 GHz band with a similar antenna configuration at a future 
favorable opposition (e.g., 16 Nov 2018, A = 1.04 AU). These 
observations would match the resolution of MIRO’s 1.6 mm 
channel during Rosetta \s flyby of Lutetia. At these scales, 
measurements of the brightness temperature will provide new 
information about the surfaces of these bodies. To develop 
accurate thermophysical models, it will be important to ob¬ 
serve them at multiple (sub)millimeter wavelengths where a 
drop in emissivity to values of ~ 0.6-0.8 has been reported 
(e.g., Muller & Lagerros 1998) particularly for rockier bodies 
(e.g., Gulkis et al. 2010). ALMA can also potentially measure 
the mutual orbit of smaller binary asteroids, providing impor¬ 
tant information on the mass and density of such objects (e.g., 
Carry et al. 2015). Finally, the ability of ALMA to deliver 
very accurate astrometry will enable better long-term model¬ 
ing of asteroid orbits, leading to improved predictions (Busch 
2009). 


This paper makes use of the following ALMA data set: 
ADS/JAO.ALMA#2011.0.00013.SV. ALMA is a partner¬ 
ship of ESO (representing its member states), NSF (USA) 
and NINS (Japan), together with NRC (Canada), NSC and 
ASIA A (Taiwan), and KASI (Republic of Korea), in cooper¬ 
ation with the Republic of Chile. The Joint ALMA Obser- 
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vatory is operated by ESO, AUI/NRAO and NAOJ. The Na¬ 
tional Radio Astronomy Observatory is a facility of the Na¬ 
tional Science Foundation operated under cooperative agree¬ 
ment by Associated Universities, Inc. This research has made 


use of NASA’s Astrophysics Data System. We thank Thomas 
Muller, Mark Gurwell, Bryan Butler, Rafael Hiriart, Ralph 
Marson, Dirk Petry, and Vivek Dhawan for useful discussions. 
Facilities: ALAiK. 
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